Geometric phases of mesoscopic spin in Bose-Einstein condensates 
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We propose a possible scheme for generating spin- J geometric phases using a coupled two-mode 
Bose-Einstein condensate (BEC). First we show how to observe the standard Berry phase using 
Raman coupling between two hyperfine states of the BEC. We find that the presence of intrinsic 
interatomic collisions creates degeneracy in energy that allows implementation of the non-Abelian 
geometric phases as well. The evolutions produced can be used to produce interference between 
different atomic species with high numbers of atoms or to fine control the difference in atoms 
between the two species. Finally, we show that errors in the standard Berry phase due to elastic 
collisions may be corrected by controlling inelastic collisions between atoms. 



I. INTRODUCTION 

Geometrical phases in quantum theory have attracted 
considerable interest since Berry [|] showed that the state 
of a system undergoing an adiabatic and cycUc evolu- 
tion acquires this purely geometric feature in addition to 
the usual dynamical phase. The term "geometric" comes 
from the fact that the phase factor acquired by the state 
depends only on the path followed by the state but not on 
the rate at which it is traversed. If the system is nonde- 
generate, the geometric phase is simply a complex num- 
ber called the Abelian phase, but in general it is a unitary 
matrix inducing transitions between degenerate states, 
which is called the non-Abelian phase or holonomy 0. 
Geometric phases have been proposed [|| and tested |0] 
in a variety of settings. Generalizations of Berry's phase 
analysis to nonadiabatic, noncyclic, and nonunitary evo- 
lutions have also been achieved [^. Sjoqvist et. al. pro- 
posed an operationally well defined generalization of the 
geometric phase for mixed states |^ and more recently, a 
fully quantized version of the phase has been given that 
considers vacuum induced effects Q. Holonomies and 
Berry phases have relevance in implementing quantum 
computation |^ , where a universal set of quantum gates 
can be performed in a fault-tolerant way by a succession 
of geometrical unitaries Moreover, geometric evolu- 
tion of states may also have importance in manipulat- 
ing quantum systems such as Bose-Einstein condensates 
(BECs). For example, a certain type of geometric phase 
have been used to create vortices in BECs jl^]. In this 
paper, we propose a method of testing both Abelian and 
non-Abelian geometric phases for a general spin- J sys- 
tem modeled by two coupled Bose-Einstein condensates. 
This is interesting, because it allows us to test geometric 
phases for spin values as mesoscopic as J ^ 10*, a feat 
not accomplished yet in any other system. 

Trapped atomic BECs, first achieved some seven years 
ago, provide us with the ability to make mesoscopic quan- 
tum objects containing of the order 10^ atoms in the same 



quantum state. Although a mesoscopic system may, in 
principle, be constructed using photons as well, a crucial 
distinction is that atoms may be stored for longer times. 
A longer storage time (of the order of an hour for mag- 
netic traps |pl]|) implies a longer decoherence time scale 
of quantum states, thereby aiding the implementation of 
adiabatic evolution (required for the Berry phase). Ge- 
ometrical phases are generated by adiabatically varying 
the Hamiltonian of a system in a cyclic fashion. This 
can be done in BECs since the states of ultracold atomic 
samples can be manipulated by electromagnetic fields. 

The Hamiltonian describing a fixed number of atoms 
in two different internal levels trapped by a magnetic po- 
tential can be approximated, for systems composed of a 
few thousand atoms, by a two- mode Hamiltonian p^ , p^ . 
The Schwinger oscillator model allows to conveniently ex- 
press the two-mode problem in terms of angular momen- 
tum operators. In this way the Abelian and non-Abelian 
phases generated are associated with spin- J states, where 
the spin is related to the total number of atoms in the 
condensate. We use two-photon excitations, which gen- 
erate a coherent superposition in the two- mode BEC, to 
demonstrate the Berry phase when the collisions between 
atoms are neglected. We then extend this idea further 
to consider the possibility of measuring holonomies in- 
cluding the nonlinear term due to coUisional interaction 
between particles for which the Hamiltonian is degen- 
erate. To produce Abelian and non-Abelian transfor- 
mations when the non-linear term in the Hamiltonian 
is considered, in addition to the two-photon excitations, 
inelastic collisions must be manipulated. Together the 
two-photon excitations and the inelastic collisions cor- 
respond to a two-mode displacement of the eigenstates 
of the system and this may have potential applications 
in manipulating BECs. We shall show that, if collisions 
between atoms are seen as errors in the Berry phase gen- 
erated by the linear Hamiltonian, these errors can be 
corrected by introducing inelastic collisions. 

The paper is organized as follows. In Sec. II we de- 
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scribe, using angular momentum operators, the physical 
system we consider for generating geometric phases: the 
two-mode BEC. Section III is devoted to explain how 
Berry phases and holonomies arise, presenting the main 
mathematical formalism used through the paper. In sec- 
tion IV we discuss how to generate Berry phase in the 
BEC described in section II and we propose a scheme 
for detecting this phase. In Sec. V we move to consider 
the holonomies related to the two-mode BEC. Section VI 
contains final remarks concerning the generation of the 
Berry phase when collisions between the atoms in the 
condensates cannot be neglected, and finally the paper 
concludes in Sec. VII. 



II. THE TWO-MODE BOSE-EINSTEIN 
CONDENSATE 

A physical realization of our system is two condensates 
in different hyperfine levels \A) and \B), such as those al- 
ready produced in experiments by Myatt et. al and 
Stenger et. al . In the Rb system of Ref. , an ex- 
ternal laser is applied to induce a Josephson-like coupling 
and the detuning of the laser is adiabatically changed 
to produce various transitions. Alternatively, for the Na 
spinor system of Ref. , state-dependent magnetic field 
gradient may be applied to induce Josephson tunneling. 
The Hamiltonian for the system can be written under the 
two-mode approximation, taking the annihilation opera- 
tors to be a and h for the two distinct hyperfine states: 



H 


— Ha + Hb + Hint + Hlas, 


(1) 


Ha 


t , t t 
= LOaO, a H — 2' 


(2) 


Hb 


= ojbb^b+^b^b^bb, 


(3) 


Hint 


— ^^^a^ ab^b, 
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(4) 


Hlas 




(5) 



where Ha and Hb describe the two condensates undergo- 
ing self-interactions and Hint and Hias describe the con- 
densates interacting with one another via coUisional and 
laser-induced interactions, respectively. A is the detun- 
ing of the laser from the |^) — > \B) transition. We note 
in particular that Hias describes Josephson-like coupling 
which interchanges internal atomic states in a coherent 
manner. 

As discussed above, the Hamiltonian can be written 
in a more suitable way by employing the Schwinger an- 
gular momentum [S'J7(2)] operators defined as Jx — 
i(at6-Ha5t), = {l/2i){a'<b-ab^),andJ, = ^{a''a-b^b). 
The Casimir invariant J'^ = + Jy + Jl has eigenval- 
ues i{i -1-1) and j represents the total number 2A^ = 
2{Na + Nb) of the two different species of atoms {N 
could be 10^ for the type of condensates we consider). In 
terms of Schwinger operators, the Hamiltonian takes the 
simple form 



H = aJ,+ PJI -I- 7[cos((/)) -I- sm{4>)Jy\ (6) 

where a = tj,, - Wf, + (2 J - l){Ua - Ub)/2, p ^ {Ua + Ub- 
Uab)/2, and (j) = At. The eigenvalues m of the opera- 
tor Jz represents the difference 2{Na — Nb) in the num- 
ber of atoms in different hyperfine levels, while Jx (and 
Jy) takes on the meaning of the relative phase between 
the two species. It is noted here that s-wave scattering 
lengths may be tuned using Feshbach resonances by the 
application of an external magnetic field The fac- 

tors Ua, Ub, Uab, and consequently a and (3 are therefore 
adjustable parameters giving us important additional de- 
grees of freedom. The above Hamiltonian [Eq.(^)] is suffi- 
cient to generate the Abelian Berry phase if the collision 
term (/JJ^) is negligible. We shall describe later how 
terms such as JxJy and (required for the non- Abelian 
geometric phase) can be added to the above Hamiltonian 
by inducing inelastic collisions. 



III. GEOMETRIC EVOLUTIONS 

We will now describe briefly how geometric phases 
arise. Holonomies are unitary transformations of geo- 
metric origin generated by varying a set of external pa- 
rameters a = {c/ilj^^o featuring in the Hamiltonian of a 
system, in a cyclic and adiabatic way. Their dimension- 
ality n equals the degree of degeneracy of the eigenspace. 
The Berry phase is the special case when the eigenspace 
is nondegenerate, and the unitary transformation is then 
one dimensional, i.e., a complex number. Consider that 
the adiabatic variation of the Hamiltonian is given by 
H{a) = U{a)HQU^ {a) (where the parameters a vary on 
some control manifold AA) with the lA{a) transformation 
being such that the degeneracy structure of the initial 
Hamiltonian i?o is preserved. For this purpose, the pa- 
rameters must be varied slowly with respect to any time 
scale associated with the system dynamics. After a com- 
pletes a loop C in A^, an initially prepared state \^in) 
is mapped to \^out) = e'^'^^'^T a{C)\^ ^n) where T is the 
overall time of the evolution, and E is the energy of the 
degenerate space where |^m) and \^out) belong. T a{C) 
(termed the holonomy associated with non- Abelian con- 
nection forms) appears due to the nontrivial topology 
structure of the degenerate space and is given by 

r^(C) =Pexp (a. (7) 
Jc 

P is the path ordering symbol and the Wilczek-Zee con- 
nection A H] is defined as 




for i,j = l,...,n parametrizing states belonging to the 
same degenerate eigenspace of Hq . 
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An alternative interpretation of the holonomy, which 
is needed here to calculate explicitly the integral in 
Eq. (0), is that the holonomy can be seen as the ex- 
ponential of the flux of the field strength F^^ {a, t) = 
—daAr+d-rA^ + lA^, At], through a surface parametrized 
by (T and r. If we consider the surface to be a rectangular 
loop C in the plane (u, r) with ordered sides 1, 2, 3 and 
4 , and represent the path order exponential integrals of 
the connection for each side by Wi for i — 1,...,4, for 
T~^{(7,t) = W4W3 we can state the non-Abelian Stokes 
theorem as 

TAiC) = W4W3W2W1 = PrfAiO, T^'P'^Tdadr ^ 

where P,- is the path ordering symbol with respect to 
the r variable only, unlike the usual path ordering sym- 
bol P, which is with respect to both variables a and r. 
The evaluation of the holonomies here is performed by 
the application of (^. 

IV. SPIN J BERRY PHASE 

To generate a Berry phase in the BEC, let us first ig- 
nore the terms in the Hamiltonian due to collisions. A 
BEC can be described by such a Hamiltonian when one 
may assume that the condensate is dilute enough that the 
coUisional interactions between the atoms (both inter- 
and intraspecies) become insignificant compared to the 
coupling rate generated by the strong external field. In 
principle, one may also consider the case where through 
the Feshbach resonance Ua + Ub = Uab {P = 0), or one 
can assume Uab = Ua = Ub, and tune the scattering 
length of either one of the species so that either Ua or 
Ub is reduced down to zero energy fl^ . Experimentally, 
one needs to ensure that losses due to three-body recom- 
binations are minimized. In this case the nonlinear term 
in the Hamiltonian vanishes so that 

H = aJz+ 7[cos(0) Ja; -I- sin{(j))Jy]. (10) 

This Hamiltonian corresponds to the motion of a spin 
J particle in a magnetic field with amplitude B — 
\/ o? +7^, whose direction is slowly varying as changes. 
= Ai can be made to vary arbitrarily slowly (so that 
the adiabatic approximation holds true) by choosing ar- 
bitrarily small detuning A. For a fixed 7, varying over 
a complete loop produces a nontrivial Berry phase. We 
now proceed to calculate this. The Hamiltonian can be 
written as 

H{<t>,e)^u{4>,e)HoU\4>,e) (ii) 

where Hq — a^Jz and IA(4>,9) = exp(— ic/ij^) exp(— i0 Jj,) 
with a — cos 9 and 7 = ao sin 9. This unitary transfor- 
mation on Hq corresponds to a two-mode displacement 
with amplitude 9/2 and phase (j). The eigenvectors of Hq 
are first rotated through an angle 9 on the xz plane and 



then rotated through an angle (f> in the xy plane (as shown 
in Fig.|l|) to obtain the eigenvectors of H. The eigenstates 
of Hq are the usual angular momentum eigenstates with 
Jz|j, to) — toIj, m). This implies that the eigenstates of 
our original Hamiltonian H are — U{4),9)\j,m) and 
we can then calculate the Berry phase using expressions 
and (^) for i = j = m. The connection components 
generated by the transformation U (0, 9) for the eigen- 
states of Hq are given by 

A^ = {j,m\UHct>,0)^U{ct>,9)\j,m) 

= -i{j,m\{sm{9)Jx + cos{9) Jz)\j , m) 

= —imcos9 (12) 

and 

Ag = {j,m\U\<j>,9)^U{<l,,9)\j,m) 

= -i{i,'m\Jy\j,'m) =0. (13) 

The field strength is F^^ — imsm9 and wc find that, 
if we vary 9 from zero to a fixed value and <j) from to 
27r, the state | j, m) acquires a phase 

<t>Berry=i / F^^g, d9' d(p ^ ~2Tim{l - COS 9) , {U) 

Jo Jo 

which is TO times the solid angle subtended by the cir- 
cuit at the origin in parameter space. During the evo- 
lution, the state also acquires a dynamical phase 72^„ = 

Jq {j, ^\Ho\ji 'm) = maoT. This phase can be eliminated 
by choosing an adequate evolution time T for which the 
dynamical phase is a multiple of 2tt. Note that the Berry 
phase does not depend on T but only on the geometry of 
the loop C. 

Now we shall discuss a scheme for detecting the Berry 
phase. We first prepare the system in the state 
which, when represented in terms of the population of the 
two modes, is simply |A^) (g) |0). We then switch on laser 
fields with detuning A = 0, and vary a and 7 slowly in 
such a way that \ j,j) (following the Hamiltonian adiabat- 
ically) evolves to Z//(0,7r/2)|j, j). Now we implement an 
adiabatic loop of the Hamiltonian in the parameter space 
given by the transformation Z^t(o, 9)U{2tt, 0)l({0, 9), with 
cos9 — 1/2. The middle transformation Z/^(27r, 0) is ob- 
tained by switching on the detuning A and letting the 
Raman transition on for a time such that t = If 
we choose the time of the loop of the Hamiltonian such 
that the dynamical phase is eliminated (i.e., a multiple of 
2tt) for all states, then the evolution of the state during 
this loop will be purely due to geometric phases. Next, 
the transformation ^^'''(0, 7r/2) is applied to the state. For 
our choice of 9, the Berry phase will be such that the fi- 
nal state after all transformations will be orthogonal to 
\j, j) . The presence of a Berry phase can now be verified 
by measuring the population of the second mode, which 
will now always be nonzero. This method is a general- 
ization of the usual Hadamard-Berry-Hadamard method 
used for the detection of the spin- 1/2 Berry phase. 
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J. 




FIG. 1. The spin-J Berry phase is obtained by rotating the 
eigenvectors of Hq first through an angle 9 in the xz plane 
and then through an angle in the xy plane. The phase is 
proportional to the area traversed by the vector. 



V. HOLONOMIC EVOLUTION IN THE 
TWO-MODE BEC 

A geometric evolution in the condensate becomes a 
more complicated transformation when collisions are con- 
sidered. The nonlinear term in the Hamiltonian allows 
for the possibility of degeneracy and, by slowly varying 
the Hamiltonian through the parameters, we can then 
generate holonomies. For performing a holonomic evo- 
lution the state of the system must be confined to a 
degenerate subspace at all times. A two-dimensional 
degenerate subspace can be created by making two of 
the eigenvalues of the Hamiltonian Hq — uqJz + (3oJ'^ 
equal. Choosing ao//?o — — (2r7i + 1), the states \j,m) 
and I j, m + 1) have the same energy. We shall trans- 
form the nonlinear Hamiltonian with the same unitary 
transformation that we used for the linear one, lA{<j), 9) = 
exp{—i(f)Jz) exp{-'i9Jy). This will result in a transformed 
Hamiltonian that includes inelastic collisions, but the 
states are transformed in the same manner as in the 
linear case. In particular, we obtain U(9)HoW (9) — 
aa[sm{9)J^ + cos{9)J,)] + Po[sin^i9)J^ + cos^{9)Jl + 
sin(6') cos{9){2JxJz + 1)]. The term JxJz includes inelas- 
tic collisions of the form a'^ba'^a, while includes terms 

of the form b^. Such inelastic coUisional terms de- 
scribe processes in which atoms exchange the hyperfine 
states during a collision. We do not have to introduce 
such terms artificially in our system as such processes are 
normally present and only deliberately suppressed when 
one tries to produce a BEC. The reason for this is that 
the atoms are generally confined in a magnetic trap and 



such spin flips result in the loss of the atoms. However, 
by using an optical confinement of the BEC (l8| such 
problems are removed, as the optical dipole force is not 
selective about the hyperfine states of the atoms. The 
measures to suppress inelastic collision can be removed 
with an optical trap and, indeed, one may even enhance 
these processes by inducing suitable Zeeman shifts. This 
is possible provided the total angular momentum and en- 
ergy are conserved on collision. By freeing the channel 
through which excess angular momentum is translated 
into an overall relative rotational motion of the colliding 
atoms, the inelastic coUisional processes can be enhanced. 
By taking 9 small it is possible to meet the experimen- 
tal values for the production ratios of those terms in a 
two- mode BEC. The connection components generated 
by the transformation U (0, 9) and for the above degen- 
erate states are given by 



■ smf 



f sine 



(m + 1) cos 6* 



A> 



_ p 

2 



1 
-1 



where p = ^ (j — in) (j + m + 1) . We now consider 
the case m ~ which corresponds to almost equal 
numbers of particles in both condensates. In this case 
P ~ \/ jU + !)• For ^ large number of atoms p 3> 1, 
we can neglect terms that are small compared to p in 
the following analysis. As the connection components 
and Ag do not commute with each other we have 
to employ the non-Abelian Stokes theorem to evalu- 
ate the holonomy Ta(C). Indeed, following the proce- 
dure presented in [ p9| , p0| , for a rectangular loop C with 
vertex coordinates {(0, ^o), (0 = '''/(psineo), ^o), ((/i = 
7r/(psin6'o), 6*1), (0, 6*1)} (shown in Fig.||) we obtain the 
holonomy 



exp 



-i (cos 6*1 — cos 00 ) 



(m 



sm f^o 



■ (T2 



where a2 is the Pauli matrix. To obtain the above re- 
sult we have used the approximation that for 9q large 
compared to the variation 9i ~ 9o, the function sin(e) 
is almost constant compared to sin{p9). As an appli- 
cation for 9i — 9q = 7r/(2p) we can obtain a change of 
state from \j, m) to \j, to -t- 1) (transfer of two atoms 
from one mode to the other). In addition, starting with 
the atom number state |iV^) (g) liV^) and for 9i ~ 9q — 
7r/(4p), one can produce up to an overall phase the state 
{\Na)<E)\Nb)-\Na + 2)(E)\Nb-2))/V2. This interference 
procedure can be used, e.g., in high-precision measure- 
ments for the construction of quantum gyroscopes pll . 
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{OA) 



(OA) 



(7r/psin 6*0, 6*1) 



(7r//9sin 6*0, 6*0) 



FIG. 2. The loop C for the non-Abehan Stokes theorem. 



VI. FINAL REMARKS 

The standard Berry phase can also be generated when 
collisions are included. By relaxing the degeneracy condi- 
tion ao/Po = — (2to+ 1) and performing the same trans- 
formation U ((/), 6) to the Hamiltonian the same Berry 
phase as in Eq. (|l^) is generated even in the presence of 
collisions. This is, however, true only when both elastic 
and inelastic collisions are considered. One can think of 
this in the following way: Considering elastic collisions to 
be errors, in order to generate the same Berry phase as in 
the collision-free Hamiltonian, there must be a correction 
achieved by including inelastic terms. 

Finally, we shall give an analytic expression for the 
unitary transformation that produces from the nonlinear 
Hamiltonian Hq the term "fj^ for small 7, but arbitrary 
ao and (3o- The only condition needed in order for the 
perturbation expansion to be valid is that the Hamilto- 
nian should not be degenerate, i.e., ao/Po should not be 
an integer. Indeed, with easy algebraic steps we have 



H^U{cf>)U{j)HoU^{j)U\cb) 



(15) 



where the dependence is exact, while the 7 dependence 
is valid for weak Josephson-like coupling, 7^1. The 
unitaries are given by U{(l>) — e:sq){i(j>Jz) and C/(7) = 
exp(i7G) with 



G ^ a,ki Jz'^ Jy Jz'' 



a-ki = (-1) 



kill a^+'+i 



U{(f>) and U{'-f) rotate the basis of eigenvectors of Hq. As 
G is evaluated only to first order in 7 it is impossible 
to evaluate from it a nonvanishing Berry phase, as its 
calculation involves two exterior derivatives of the trans- 
formation U{'y). 



VII. CONCLUSIONS 

We have presented a procedure for evolving the state 
of a two-mode BEC in a geometrical fashion. The BEC 
setup described here due to its two mode bosonic na- 
ture resembles other bosonic models studied in connec- 
tion with Berry phases and holonomies Here we 
have focused on the spin-J description of the conden- 
sate that is apparent during the measuring procedure 
of the Berry phase. In particular, in the limit where 
collisions can be neglected the state of the condensate 
acquires a Berry phase by varying the displacement pa- 
rameter in a cyclic and adiabatic way. This resembles 
a spin-J particle, where J is mesoscopic. Berry's phase 
is manifested by varying the direction of the magnetic 
field. Considering the degeneracy introduced by the col- 
lisions between atoms, a holonomic evolution, generated 
by slowly changing the parameters, allows for controlled 
transfer of population between modes. In addition to al- 
lowing for tests of Abelian and non-Abelian geometrical 
phases for mesoscopic J, this might also be useful as a 
procedure for manipulating the condensate. 
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